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ABSTRACT 

Using long-slit optical spectroscopy obtained at the 10.4 m Gran Telescopio Canarias, 
we have examined the gaseous environment of the radio-loud quasar TXS 1436+157 
(— j ■ (z=2.54), previously known to be associated with a large Lya nebula and a spatially 

O [ extended Lya-absorbing structure. From the Lya nebula we measure kinematic prop- 

erties consistent with infall at a rate of ^10-100 Af© yr _1 - more than sufficient to 
power a quasar at the top of the luminosity function. 

The absorbing structure lies outside of the Lya nebula, at a radius of >40 kpc 
from the quasar. Against the bright unresolved continuum and line emission from 
the quasar, we detect in absorption the NV AA1239,1241, CIVAA1548.1551 and SilV 
AA1394,1403 doublets, with no unambiguous detection of absorption lines from any 
low-ionization species of metal. The metal column densities, taken together with the 
HI column density measurement from the literature, indicate that the absorbing gas 
is predominantly ionized by the quasar, has a mass of hydrogen of > 1.6 x 10 11 Mq, 
a gas density of ^18 cm -3 , a line of sight thickness of >18 pc, and a covering factor 
approaching unity. While this absorbing structure is clearly not composed of pristine 
gas, it has an extremely low metallicity, with ionization models providing a 3er limit 
of 12+log(0/H)<7.3. 

To explain these results, we discuss a scenario involving starburst-driven super- 
bubbles and the creation of infalling filaments of cold gas which fuel/trigger the quasar. 
We also discuss the possibility of detecting large-scale absorbers such as this in emis- 
sion when illuminated by a powerful quasar. 

Key words: quasars: emission lines; quasars: absorpion lines; quasars: individual: 
TXS 1436+157; galaxies: ISM; galaxies: active; galaxies: evolution; ultraviolet: galax- 
ies; ultraviolet: ISM 



1 INTRODUCTION 

Lya nebulae at high redshift (z>2: HzLAN), sometimes also 
known as Lya 'blobs', 'halos' or 'fuzz', are prodigious sources 
of HI Lya photons, with luminosities up to ~10 45 erg s _1 , 
and sizes often exceeding ~100 kpc (e.g. Steidel et al. 2000). 
HzLANs are not associated with a common type of galaxy: 
they are variously associated with UV-selected star forming 
galaxies, dusty star-forming galaxies, powerful active galac- 
tic nuclei, or with no apparent optical galaxy at all (e.g. Mc- 
Carthy et al. 1987; Reuland et al. 2003; Webb et al. 2009; 
Steidel et al. 2011; Bridge et al. 2012). A complete under- 
standing of the nature of HzLANs has not yet been reached, 
with key properties such as their dominant power source(s) 
and their space density remaining particularly poorly un- 



derstood. Nevertheless, HzLANs are thought to represent an 
important phase of mass assembly, and understanding this 
energetic phenomenon promises to yield significant insights 
into the physics of massive galaxy formation (e.g., Mori & 
Umemura et al. 2006) . 

Associated absorbers offer a complementary means by 
which the gaseous environment of galaxies can be stud- 
ied (e.g. Pettini et al. 2002; Fosbury et al. 2003; Wilman 
et al. 2005). A substantial fraction of radio-loud, powerful 
active galaxies at high-z show spatially extended, narrow 
(FWHM<500 km s" 1 ) absorption lines, usually HI Lya, 
close to the redshift of the galaxy; the observational char- 
acteristics of the absorbers suggest they are giant shells of 
gas, with radii >50 kpc, and are distinct structures from 
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the narrow line emitting nebulae that are also usually asso- 
ciated with powerful active galaxies (Rottgering et al. 1995a; 
van Ojik et al. 1997; Binette et al. 2000; Jarvis et al. 2003; 
Wilman et al. 2004; Humphrey et al. 2008a). Intriguingly, 
in the case of high-z radio galaxies there is a strong anti- 
correlation between the projected size of the radio source 
and the detection of associated HI absorbers, with ~90 per 
cent of high-z radio galaxies that have small radio sources 
(<50 kpc) showing an associated absorber, falling to ~25 
percent for large radio sources (>50 kpc: van Ojik et al. 
1997). This result has been interpreted in terms of (a) cre- 
ation of an expanding shell of gas swept up from the inter- 
stellar medium of the host galaxy by the expanding radio 
source (e.g. Krause 2002), or (b) pre-existing gaseous shells 
produced by powerful starbursts, which are over-taken and 
disrupted by the expanding radio source (e.g. Jarvis et al. 
2003). It is not yet clear whether the absorbing shells are 
in infall or outflow (e.g. Jarvis et al. 2003; Humphrey et al. 
2008a) . At least a few of the HI absorbers also contain car- 
bon, as shown by their detection in CIV A1549 absorption 
lines (e.g. Binette et al. 2000; Jarvis et al. 2003; Humphrey 
et al. 2008a). Binette et al. (2000, 2006) and Jarvis et al. 
(2003) have used ionization models to place important con- 
straints on the metallicity, the source of ionization and the 
density of the HI and CIV absorbers detected in front of 
radio galaxies USS 0200+015 (z=2.23) and USS 0943-242 
(z=2.92). Interestingly, Vernet et al. (2001) detected vari- 
ous narrow absorption lines from low-ionization or neutral 
species of C, O, S and Si against the UV continuum emis- 
sion of several high-z radio galaxies, and noted a tendency 
for the absorption lines to be stronger in the galaxies with 
more highly polarized UV continuum emission. The possible 
relationship between these low ionization metal absorption 
lines and the spatially extended HI absorbers is not yet clear. 

Quasars provide a useful alternative perspective from 
which to examine narrow associated absorbers. The rela- 
tively unobscured view to their highly luminous central en- 
gines significantly enhance the detectability of narrow ab- 
sorption lines that would, as in type II QSOs or radio galax- 
ies, otherwise lie in spectral regions with no strong contin- 
uum or line emission. The absorbers also ought to have a rel- 
atively unobscured view of the AGN, raising the possibility 
that some of the narrow absorbers associated with quasars 
may be photoionized by the hard radiation field of the AGN; 
this may result in the detection of absorption lines from more 
numerous metal species than neutral or stellar-photoionized 
absorbers, greatly aiding analyses of excitation and metal- 
licity. Moreover, by examining the properties of associated 
absorbers as a function of AGN and jet orientation, informa- 
tion about the geometry of the absorbers can be obtained. 
Indeed, Baker et al. (2002), in their study of associated ab- 
sorption in radio-loud quasars, found that the detection of 
associated narrow CIV absorbers is orientation-dependent, 
being detected more frequently in the lobe-dominated or 
steep spectrum quasars than in the flat-spectrum or core- 
dominated quasars, suggesting that the absorbing gas lies 
away from the axis of the radio jet. The authors also noted 
an anticorrelation between the equivalent width of narrow 
CIV absorption lines and the projected diameter of the radio 
source, broadly similar to the trend identified by van Ojik 
et al. (1997) for high-z radio galaxies. 

We have begun a programme of spectroscopic observa- 



tions, designed to help elucidate the nature of Lya nebulae 
and HI absorbers associated with powerful active galaxies at 
high-z. By targeting active galaxies known to show extended 
Lya nebulae, extended HI absorption close to the redshift of 
the galaxy, and a very bright unresolved emission component 
(>10 -17 erg s cm~ 2 A" 1 ), we aim to obtain the geometric, 
ionization and chemical properties of individual absorbers. 
This approach also mitigates some problems inherent to pre- 
vious studies, such as the lack of sufficiently bright back- 
ground emission to readily detect absorption lines other than 
HI Lya or CIV, or the lack of an extended background source 
to provide information about the size and location of the 
absorber. Based on this idea, Humphrey (2008a) employed 
deep longslit spectroscopy from Keck II, together with in- 
tegral field spectroscopy from the Very Large Telescope, to 
examine the properties of a spatially extended HI absorber 
associated with the reddened radio-loud quasar MRC 2025- 
218. The authors traced the HI absorption feature across a 
projected area of ~30 kpc x 40 kpc, with very little varia- 
tion in kinematic properties therein, which implies the ab- 
sorber is a coherent structure located at least 20 kpc from 
the AGN. Against the unresolved continuum and line emis- 
sion from the active nucleus, narrow absorption lines of CI, 
CII, CIV, NV, OI, Sill, SilV, A1II and A1III were detected 
at the redshift of the HI absorber. The relative column den- 
sities of the metal absorption lines indicate the absorbing 
gas is photoionized by the active nucleus, has a gas density 
of >10 cm -3 , and is enriched in metals. 

Following on from Humphrey et al. (2008a), in this pa- 
per we present long-slit spectroscopic observations obtained 
at the 10.4 m Gran Telescopio Canarias of the R — 19.6 
quasar TXS 1436+157 at z=2.537. With an 4.7 GHz flux 
density of 47 mjy, this quasar qualifies as radio-loud, and its 
radio source shows a steep spectrum and a lobe-dominated 
morphology (Rottgering 1993; Carilli et al. 1997). Of par- 
ticular relevance to this programme, TXS 1436+157 is also 
known to be embedded within a ~90 kpc diameter nebula 
with a Lya luminosity of ^ 2.2xl0 44 erg s _1 (van Ojik et 
al. 1997: Rottgering et al. 1997), and shows a strong and 
spatially extended Lya absorption feature within its Lya 
emission line profile (van Ojik et al. 1997). The primary 
aim of this paper is to elucidate the properties of the Lya 
nebula and the associated absorber. We assume Ho=71 km 
s" 1 Mpc~\ fi A =0.73 and fi m =0.27. At the redshift of TXS 
1436+157, l"corresponds to 8.2 kpc. 



2 OBSERVATIONS AND DATA REDUCTION 

Long-slit spectroscopic observations of TXS 1436+157 were 
obtained in service mode on 2010 May 16 and May 21, us- 
ing the Optical System for Imaging and low Resolution In- 
tegrated Spectroscopy (OSIRIS) instrument at 10.4 m Gran 
Telescopio Canarias. The observing conditions were dark 
and photometric. The total integration of 7590 s was split 
into 3 integrations of 2530 s each, one of which was car- 
ried out on 2010 May 16, while the other two were carried 
out on 2010 May 21. Although the FWHM of the seeing 
disc was l"during all of the observations, the primary mir- 
ror segments were not optimally aligned during the observa- 
tions, which resulted in a poor spatial point-spread function 
along the slit (FWHM=2.3", 1.6"and 2.3"). The R1000B 
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Figure 1. Image of the 40" X 40"field around TXS 1436+157, 
taken through a Sloan r filter with an exposure time of 20 s. The 
position of the long slit is shown. Four discrete sources were within 
the slit: the z=2.54 quasar TXS 1436+157 (Q); two galaxies which 
we know now to be at z=0.160 and z=0.159, respectively (G); and 
a continuum only source of unknown redshift (C). 



grism was used, yielding a spectral range of 3700-7000 A. A 
1.23"wide slit was used, and was oriented at a position angle 
of 50° North through East, so as to include TXS 1436+157 
and three additional sources that are within ~30"of the 
quasar, and which are closely aligned with the major axis 
of the quasar's radio emission (Rottgering et al. 1995). Fig. 
1 shows our 20 s acquisition image taken through a Sloan 
r filter, with the slit position overlaid. The instrument was 
read-out with 2x2 binning, resulting in a spatial scale of 
0.25"pixel _1 , and a wavelength scale of 2.22 A pixel -1 The 
data were reduced in the standard way using long-slit reduc- 
tion software from the IRAF software suite. Flux calibration 
was performed using observations of GRW+708247. 

Since the slit width (1.23") was significantly larger than 
the seeing disc FWHM the instrumental profile (IP) 

differs between unresolved and extended sources. For ex- 
tended sources the FWHM of the IP is 11.0 A, measured 
from night sky lines in the spectra; this IP will be adopted 
when calculating the FWHM of emission lines from spatially 
extended line-emitting gas. For unresolved sources, it is ex- 
pected to be ~8 A, in good agreement with narrow emission 
lines in the spectra of the intervening star forming galaxies 
and from the strong absorption lines in the spectrum of the 
quasar; this IP will be used for absorption lines detected 
against the unresolved line and continuum emission of the 
quasar. 



3 INTERLOPING SOURCES 

The broad band image of TXS 1436+157 presented by 
Rottgering et al. (1995b) shows four tightly aligned sources, 
including the active galaxy, running along a position angle 
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Figure 2. Sections of the one-dimensional spectrum showing the 
absorption lines detected in front of TXS 1436+157. In the lower 
panel, the gap in the data between CIV and Hell hides a sub- 
stantial sky-subtraction residual. 
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Figure 3. Sections of the one-dimensional spectra of TXS 
1436+157 and the two intervening HII galaxies, with the ex- 
pected wavelength of the OI A1302.2 absorption line marked with 
a dashed line. The flux scale is for TXS 1436+157; the spectra of 
the intervening galaxies have been scaled up and smoothed. Note 
that all three of the galaxies show a strong absorption feature 
close to this wavelength. 

of ~40° (see our Fig. 1). The radio source shows a roughly 
similar position angle (Carilli et al. 1997). Our long slit was 
positioned so as to encompass all four of these sources, in or- 
der to investigate the nature of these aligned optical sources. 
Our long slit spectrum reveals that the two closest sources 
on the sky to TXS 1436+157 are in fact emission line galax- 
ies at intermediate redshifts. At a distance on the sky of 
3" NW, SDSS J143904.77+153121.5 has z=0.159, while 
SDSS J143905.47+153112.6 at a distance of 11" SE has 
z=0.160. Thus, the alignment of these nearby (in projec- 
tion) sources with the radio source is merely coincidental. 
In the case of the source 25" to the SE, we do not detect 
any strong absorption or emission features, and thus we can- 
not ascertain its redshift. 



4 RESULTS AND DISCUSSION 

4.1 Measuring W\ for absorption lines 

For a formal detection we require (i) that the FWHM of 
the candidate absorption feature is not smaller than that 
of the instrumental profile, (ii) that difference between the 
continuum level and the minimum of the feature is at least 
twice the local la pixel to pixel uncertainty in flux density, 
and (iii) in the case of doublets, that the wavelength separa- 
tion and W\ are in agreement with their theoretical values. 
Since they are doublets and satisfy criterion (iii), our iden- 
tifications of the NV, SilV and CIV lines are rather secure 
(>10er). Due to the low spectral resolution of the data, ab- 
sorption lines were measured using single negative Gaussian 
profiles, in order to obtain their equivalent width {W\). Un- 
certainties are derived from uncertainty in the level of the 
underlying continuum and broad line emission. 

The CIV absorption doublet has altered the overall ve- 



Table 1. Measurements and limits on the absorption lines. 
Columns: (1) species; (2) rest-frame wavelength (A); (3) velocity 
shift (km s~ 1 ) relative to the central wavelength of the nuclear 
narrow Hell A1640 emission line; negative values indicate a blue- 
ward shift; the la uncertainties are 50-70 km s ; (4) observer- 
frame equivalent width (A); we show for completeness the lower 
limits for HI Lya and CIV; (5) oscillator strength of the transi- 
tion^); (6) column density of the atom/ion. The rows beginning 
'NV mean' and 'SilV mean' give 1/a 2 weighted means for iVjv_|_4 
and Ngi+s, respectively. 'Sill total' gives a combined upper limit 
for the 1260.4 A line and the 1264.7 A excited fine structure 
line. ^The HI velocity and column density is taken from van Ojik 
ct al. (1997). 
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locity profile of the CIV line to such an extent that the veloc- 
ity profile of the unabsorbed emission, and hence W\, cannot 
reliably be reconstruced - at least at this spectral resolution 
(Fig. [2} . Our determine a lower limit on W\ was measured 
assuming the two emission peaks of the CIV profile repre- 
sent the maximum height of the intrinsic (unabsorbed) line 
profile. The velocity profile of HI Lya presents a similar chal- 
lenge at our spectral resolution. In this case, we refer to the 
moderate resolution spectrum and fitting presented by van 
Ojik et al. (1997), who favoured N H i ~4.0xl0 19 cm" 2 . 

Table[TJlists the detected absorption lines. Velocities are 
given relative to the centroid of the narrow Hell A1640 emis- 
sion at the position of the quasar. Also listed are upper limits 
for those absorption lines that are undetected in our data, 
but which are interesting in the context of the goals of our 
study. Fig.[2]shows regions of the extracted one-dimensional 
spectrum around the detected NV, SilV and CIV absorption 
lines. Table [TJ also includes upper limits for several poten- 
tially useful absorption lines. These limits were determined 
by calculating the equivalent width that would correspond 
to a gaussian absorption feature that has an FWHM equal 
to that of the appropriate instrumental profile (8 A), and a 
minimum that is below the continuum level by at least twice 
the local la pixel to pixel uncertainty in flux density. The 
approximate expected wavelength of the OI A1302.2 line co- 
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incides with a detected absorption line at ~4606 A, but the 
presence of an absorption feature at this wavelength in the 
spectra of the two intermediate redshift star forming galax- 
ies (Fig. [3]) suggests it may instead be due to an absorption 
system in front of all three galaxies (i.e., at z<0.16). There- 
fore, we adopt the equivalent width of this aborption line as 
an upper limit to that of 01 A1302.2. In Fig [4] we illustrate 
the conspicuous absence of various other low-ionization ab- 
sorption lines. 

4.2 Column Densities using the CoG 

We use the curve of growth (CoG) method to derive column 
densities from the observed values of W\ (Spitzer 1978). 
First, we must ascertain on which part of the CoG are the 
absorption lines. For the two clearly resolved doublets NV 
AA1239,1241 and SilV AA1394,1403 the ratio between the 
equivalent width of the blue and red line, Wuue/W re d can 
be used. Wu U e/W re d=2.0, 1.1, or -\/2, for the linear, flat, or 
damping parts of the CoG, respectively. For the NV doublet, 
the Wbiue/Wred~'2.1±0.5 ratio we measure places it on the 
linear part of the CoG. In the case of the SilV doublet the 
measured Wuue/W re d~1.6i0.6 is ambiguous. However, if 
we assume that the SilV lines have a similar intrinsic width 
to the NV lines, then the fact that Wsnv ~ Wnv would 
place the SilV lines on the linear part of the CoG also. 

For an absorption line on the linear part of the CoG, 
the column density N can be calculated using 

^ 1.13X10^,0 
/ Ao J 

where f is the oscillator strength, Wa,o is the rest-frame 
equivalent width (in A) and Ao is the rest-frame wavelength 
of the line (in A), and N is in units of cm -2 . The resulting 
a _1 -weighted column densities are N^v =2.45±0.23 xlO 14 
cm" 2 and A S i/v=0. 99+0.09 xlO 14 cm" 2 . For CIV we apply 
this relation to our lower limit to Wciv to obtain a lower 
limit N C iv ^3.8 xlO 14 cm -2 . 

4.3 Properties of the Absorber 

4-3.1 Metallicity and ionization 

What can be said about the absorber independently of pho- 
toioniation models? The detection of metal absorption lines 
from ionized species clearly demonstrates that the absorber 
is not composed of pristine material. Moveover, it must be 
at least partially ionized, with the detection in absorption of 
the NV doublet signalling photoionization by a hard radia- 
tion field with a significant photon luminosity at hv ^ 77.5 
eV (A < 160 A). Given that N C iv / Noi + Nan ^3.0 and 
Nsnv I Nsn + Nsui ^2.7, we can say that the absorber 
must be predominantly ionized. 

In order to further investigate the ionization and chem- 
ical composition of the absorber, we have computed pho- 
toionization models using the code MAPPINGS Ic (Binette 
et al. 1985; Ferruit et al. 1997). We adopt an ionizing spec- 
tral energy distribution appropriate for a powerful quasar, 
characterised by a power-law with spectral index a = —1.5 
(e.g. Robinson et al. 1987; this is justified by the detec- 
tion of NV absorption lines). The hydrogen density is set 
to uh = 10 cm -3 , but in the low-density regime, using a 



different value does not significantly affect the model re- 
sults. We adopt an isochoric behaviour for the gas density. 
We take the ionization parameteiQ U to vary from 0.0001 
to 1.6 in steps of a factor of 2. We adopt the solar abun- 
dances of Anders & Grevesse (1989), except in the case 
of oxygen for which we adopt the more recent determina- 
tion of 0/11=4.9 x 10~ 4 made by Allende Prieto, Lambert & 
Asplund (2001). We compute models for four different gas 
metallicitiefB Z=8.69 (solar), Z=7.69 (1/10 solar), Z=6.69 
(1/100 solar) and Z=5.69 (1/1000 solar). In the interest of 
simplicity, we scale all metals equally with oxygen. Each 
model calculation is ended once the HI column density has 
reached 4xl0 19 cm -2 , to be consistent with the observa- 
tions. This means at low-U the models have a large neutral 
zone, while at the very high-U end the models become mat- 
ter bounded. 

van Ojik et al. (1997) found that the absorbing gas is 
kinematically very quiescent, with a Doppler parameter of 
37±20 km s _1 measured from the Lya absorption feature. 
We have not considered ionization by shocks on account of 
this kinematic quiescence. 

Figure [5] shows various column density ratios resulting 
from our photoionization model calculations. We have plot- 
ted various column density ratios involving high-ionization, 
low-ionization and neutral species in order to break the 
degeneracy between metallicity and ionization state. Our 
lower limits on CIV/CII and SilV/Sill require U>0.01 and 
Nh >10 21 cm~ 2 . In addition, the position of the absorber 
relative to the model loci requires extremely low abundances 
of the relevant metals. More specifically, 12 + log[Si/H] 
and 12 + log[N/H] are <6.7 (<1 per cent of solar), 12 + 
log[C/H] <6.2 (< 0.3 per cent of solar). 

The OI/HI ratio gives us the slightly less sensitive limit 
12 + log[0/H] <7.3 (<4 per cent solar). It is interesting to 
note this is comparable to (or below) the oxygen abundance 
in the neutral phase in the halos of local star forming dwarf 
galaxies (e.g. Thuan et al. 2005), and also in the warm ion- 
ized phase of extremely metal-poor dwarf galaxies in the 
nearby Universe (e.g. Izotov et al. 1999; Papaderos et al. 
2006; Izotov et al. 2006; Izotov et al. 2009). 



4-3.2 Geometry and Mass 

In their moderate-resolution long-slit spectrum of TXS 
1436+157, van Ojik et al. (1997) detected the HI Lya ab- 
sorption feature out to a projected radius of ~40 kpc from 
the quasar. This provides us with a rough lower limit of 
>40 kpc for the distance of the absorbing structure from 
the quasar. If the absorbing gas is arranged in a spherical 
shell with radius r^40 kpc encircling TXS 1436+157, then 
its hydrogen mass can be estimated using 

M H = A-Kr 2 m H N H (2) 



1 We define ionization parameter as U = -. — ^ , where Q is 

the luminosity of ionizing photons of the ionizing source, r is the 
distance of the absorber from the ionizing source, and Nh is the 
hydrogen density at the front face of the gas cloud. 

2 We define metallicity as Z=12+log[0/H]. 
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Figure 4. Sections of the one-dimensional spectrum of TXS 1436+157 showing the expected wavelengths of various undetected low- 
ionization lines. The flux scale is arbitrary, and no smoothing has been performed. Upper limits on W\ are given in Table [T] 





Figure 5. Photoionization model predictions for gas photoionizcd by an active galactic nucleus. Each locus represents a sequence in U 
at four different metallicities (in the 12+log(0/H) system): 8.69 (solid line: solar value), 7.69 (dashed line), 6.69 (dash-dot line) and 5.59 
(dotted line). U increases from the left. Arrows show the measurement or limits on the relevant column density ratios. 
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Table 2. Emission line ratios measured in the extended line emit- 
ting nebula 2.0"- 3.8"East of the quasar. For comparison, we also 
show the values of the ratios from the z>2 radio galaxy (HzRG) 
average emission line spectrum of Humphrey et al. (2008b) . 



Ratio 


Value 


HzRG average 


Lya 1 Hell A1640 


22±4 


11.7 


NV A1241 / Hell A1640 


^0.3 


0.3 


NV A1241 / CIV A1549 


^0.2 


0.5 


NV A1241 / Lya 


^0.014 


0.04 


CIV A1549 / Hell A1640 


1.3±0.2 


1.8 


CIV A1549 / CHI] A1908 


4.2±1.5 


2.4 


CIII] A1908 / Hell A1640 


0.3±0.1 


0.73 



In this case the total mass of the absorber would be 
>1.6xlQ l:L Af s . If we assume instead that the absorber is a 
uniform sheet of gas then its mass would be 



Ml 



AuihNh 



(3) 



where A is the total area of the absorber. Assuming the ab- 
sorber is as extended perpendicular to the radio axis as it is 
observed to be along the radio axis (~80 kpc: van Ojik et al. 
1997), then its total hydrogen mass would be >5.1xlO lo M . 

Since we have a lower limit on U (>0.01) and on the dis- 
tance r between the quasar and the absorber (^40 kpc), and 
assuming the ionizing photon luminosity of the central AGN 
is ^ 10 57 s" 1 (e.g. Villar-Martfn et al. 2002), the density of 
the absorbing gas can then be constrained using 



n " = T^TTi ( 4 ) 

Under these assumptions, we obtain an upper limit of 
Tin ^ 18 cm -3 . The physical thickness of the absorber along 
our line of sight toward the quasar is given by 



(5) 



which gives a lower limit of I ^5.6xl0 19 cm, or ^18 pc, 
assuming a volume filling factor of 1. Significantly lower fill- 
ing factors (originally suggested by van Ojik et al. 1997) do 
not seem plausible when one considers that the absorber has 
a covering factor close to unity over several tens of square 
kiloparsecs. 

It is also interesting to estimate the mass of metals con- 
tained within the absorber. In our U>0.1 models CIV, NV 
and SilV all have an ionization fraction on the order of ~0.1, 
meaning that the total columns of C, N and Si are probably 
>3.8, ~2.5 and ~1.0 x 10 14 cm -2 , respectively. Using equa- 
tion [2l we then obtain C, N and Si masses of >6, ~4 and 
~2 xlO Mq, respectively. If we adopt the solar Si/O abun- 
dance ratio, the mass of oxygen would then be ~2 x 10 6 Mq. 
According to the calculations of Molla and Terlevich (2012), 
the production and injection of this mass of oxygen would 
require the formation of ~10 9 Mq of stars. 



4.4 Detectability of absorbing shells in Lya 
emission 

Would an absorber of the kind associated with TXS 
1436+157 be detectable in emission? For a quasar ionizing 
luminosity Q~ 10 57 s" 1 , a torus opening angle of ~ 90°, 



shell radius ~40 kpc and covering factor ~1, then the Lya 
luminosity of the absorbing shell would be ~3x 10 45 erg s" 1 , 
comparable to the Lya luminosity of the extended nebulae 
photoionized by quasars and radio galaxies. This assumes 
intervening gas in the host galaxy has a low covering factor 
and that the shell is thick enough to absorb all incident ion- 
izing photons. In narrow band imaging, such a shell would 
appear as a pair of diametrically opposed arcs centred on the 
quasar, each with a Lya surface brightness on the order of 
~10" 15 erg s" 1 cm" 2 arcsec" 2 . Thus, under the right con- 
ditions - a luminous enough quasar surrounded by a shell of 
sufficient thickness and covering factor - a gaseous structure 
of this kind should indeed be detectable in Lya emission. 

Interestingly, van Ojik et al. (1997) found that the Lya 
absorption feature of TXS 1436+157 is not black at its 
trough, with a residual flux of ~10~ 16 erg s" 1 in that spec- 
tral region. This result is qualitatively consistent with the 
idea of some in situ production of Lya photons in the ab- 
sorber, due to its photoionization by the active nucleus, al- 
though other possible explanations exist, such as a covering 
factor of <1. 



4.5 Comparison against other extended metal line 
absorbers 

In this section we compare and contrast the properties of the 
absorber associated with TXS 1436+157 the small number 
of other spatially extended metal line absorbers known to 
be associated with high-z active galaxies. 

A similar analysis to that we have performed has been 
done for another radio loud quasar, MRC 2025-218 at z=2.63 
which, like TXS 1436+157, shows a spatially extended Lya 
absorber seen in front of the galaxy's quasar-photoionized 
LAN (Humphrey et al. 2008a). In Table [3] we reproduce 
from Humphrey et al. (2008a) the properties of that ab- 
sorber. The absorbers of MRC 2025-218 and TXS 1436+157 
show remarkably similar NV column densities, which are 
3.6±1.4xl0 14 and 2.5±0.2xl0 14 cm" 2 , respectively, which 
itself suggests the two absorber systems have surprisingly 
similar (total) nitrogen column densities, and a fairly similar 

U = -. — *r ■ Moreover, both absorbers are predominantly 

ionized as shown by their large ~Nsuv /Nsir and Nc/v/Nc/ 



ratios. The limits on their gas density, although pulling in 
opposite directions, nonetheless permit a similar density of 
~10 cm" 3 . 

A comparison with the spatially extended metal line ab- 
sorbers associated with high-z radio galaxies (z>2: HzRGs) 
is also interesting. While many HzRGs show spatially ex- 
tended HI absorption in front of their extended LAN (e.g. 
van Ojik et al. 1997), only two have also been unambigu- 
ously detected in metal lines. In the case of the z=2.92 
HzRG MRC 0943-242, Binette et al. (2000) have detected 
CIV absorption at the redshift of the extended HI absorber, 
while Jarvis et al. (2003) detected in CIV absorption the 
extended HI absorber of the z=2.23 HzRG MRC 0200+015. 
Interestingly, the two HzRGs have almost identical values 
of N C iV- 3.8+0.4X10 14 and 4.9±0.2xl0 14 cm" 2 for MRC 
0943-242 and MRC 0200+105, respectively, which suggests 
their absorbers have similar (total) carbon column densi- 
ties and similar U = -. — S . Compared to these values, 

the extended absorber in front of MRC 2025-218 has a sub- 
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Figure 6. Sections of the two-dimensional spectrum of TXS 1436+157, centred on the brightest narrow emission lines: Lya, CIV, Hell 
and CIII]. The vertical scale is in arc seconds, while the horizontal scale gives observed-frame wavelength, in Angstroms. Cut levels are 
chosen to show the spatially extended, narrow emission lines, and differ between panels. The data have been smoothed by a 3 pixel 
box-car average to emphasize the extended line emission. 



Table 3. Properties of the associated absorber in front of 
the z=2.63 radio loud quasar MRC 2025-218, reproduced from 
Humphrey et al. (2008a). Av is the velocity shift of the absorber 
from the velocity centroid of the narrow Hell A1640 emission line 
measured at the projected position of the quasar. 



Species N (10 14 cm" 2 ) 

HI >300 

CI 1.4±0.3 

CII ^9.9 

CIV 20± 2 ° 

NV 3.6±1.4 

OI ^6.8 

Sil ^0.37 

Sill 13±11 

Sill* 1.7±0.7 

SilV >2.6 

A1II ^0.16 

A1III 0.4+0.1 



Property 


Value 


r 


>25 kpc 


Mh 


> 10 s M 


us 


>10 cm -3 


Av 


-600 km s" 1 



stantially higher CIV column density of N C iv=20llo x 10 14 
cm~ 2 , while the absorber in front of TXS 1436+157 has 
Nc/v ^3.8xl0 14 cm -2 which is about as high as, or higher 
than the two HzRGs. This relative enhancement of highly- 
ionized species in the extended absorbers associated with the 
quasars is consistent with the idea that those absorbers are 
photoionized by a more intense radiation field - presumably 
that of the quasar. 



4.6 The Lya nebula 

The properties of the extended Lya emitting nebula (LAN) 
are qualitatively consistent with the results presented by 
van Ojik et al. (1997). The extended Lya emission has 
FWHM=820±50 km s _1 , shows an excess of flux in its blue 
wing, and is red-shifted by 500 km s _1 relative to the Hell 
emission at the position of the quasar: these properties are 
consistent with effects expected due to the strong, spatially 
extended Lya absorption on the blue side of the line profile. 
Our new results are as follows. 



4-6.1 Kinematic properties: infall 

Out of the UV emission lines detected in the spectrum of 
TXS 1436+157, the Hell A1640 line is expected to give the 
most reliable indication of the kinematic properties of the 
gas in the LAN. As a non-resonant recombination line from 
helium its emissivity should be the least sensitive overall to 
absorption, geometry, metallicity or ionization parameter. In 
the Eastern region of the LAN the Hell line shows a signifi- 
cant red-shift of 110 km s _1 from the position of the quasar 
to >l"to the Eas10, and has FWHM=460±40 km s" 1 . These 
properties are consistent with gravitational motion within a 
massive galaxy and its dark matter halo. The quiescence of 
the kinematics are surprising given the small size of the radio 
source (45 kpc: Carilli et al. 1997), which would usually be 
associated with extreme gas kinematics (FWHM>1000 km 
s _1 : e.g., Best, Rottgering & Longair 2000; Humphrey et al. 
2006; Nesvadba et al. 2006). In order to better understand 
the kinematic properties of the Hell line, we must take into 
consideration orientation effects. 



3 On the West side of the quasar Hell emission was not detected, 
although this may simply be due to obscuration by the intervening 
galaxy at z=0.159. 
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In the standard paradigm for powerful active galax- 
ies, the ionizing radiation of the quasar is collimated into 
two diametrically opposed ionization cones by an optically 
thick toroidal structure around the nucleus. In the case of 
TXS 1436+157, our relatively direct view of the quasar then 
shows we are looking through the base of one of the ioniza- 
tion cones. Gas infalling into the host galaxy is caught in the 
ionizing radiation field of the quasar, and reprocesses ion- 
izing photons into narrow emission lines. Along the line of 
sight to the quasar, relatively bright line emission from ion- 
ized gas which is relatively closer to the quasar, and which 
is nearer to the systemic velocity of the galaxy, should dom- 
inate over low-surface brightness line emission from infalling 
gas. However, along lines of sight substantially offset from 
the position of the quasar, the line emission would be domi- 
nated by the infalling gas and will show a relative red-shift. 
If we assume that the Hell line at the position of the quasar 
is emitted close to the systemic velocity, then the relative 
red-shift of the line in the Eastern LAN indicates that the 
nebula is in inflow/infall into the host galaxy. This is con- 
sistent with previous studies which have concluded that the 
Lya nebulae around high-z radio galaxies and quasars often 
have a significant infalling/inn owing component (Weidinger, 
M0ller & Fynbo 2004; Humphrey et al. 2007; Villar-Martm 
et al. 2007a). 

Only a very rough estimate of the gas infall rate can be 
obtained from the available information. If we adopt a gas 
density of ~50 cm -2 , a filling factor ~10 -5 (e.g. McCarthy 
1993; Villar-Martm et al. 2003), and an infall velocity of 
110 km s _1 as implied by the velocity red-shift of the Hell 
line from the position of the quasar to the East region of the 
nebula, then the infall rate would be ~10-100 A/q yr" 1 . This 
is comparable to the gas infall rates estimated by Humphrey 
et al. (2007) for the giant Lya nebulae of a sample of radio 
galaxies at z>2, and is sufficient to power a quasar at the 
top of the luminosity function (Lf, a i ~10 46 erg s _1 ). 

4-6.2 Excitation: Lya excess 

To measure representative line ratios for the LAN, from 
our coadded two dimentional spectrum (shown in Fig [6]) we 
have collapsed spatially a one-dimensional spectrum from 
regions at projected distances 2.0"- 3.8"East from the po- 
sition of the quasar. The measurements were done by inte- 
grating the total flux above the continuum level, with la 
uncertainties dominated by the uncertainty in the contin- 
uum level. The line ratios are listed in Table [2] The three 
line ratios formed from only CIV A1549, Hell A1640 and/or 
CIII] A1908 are fairly consistent with the z>2 radio galaxy 
average (Humphrey et al. 2008b; see also McCarthy 1993). 
The detection of the Hell line shows that the LAN is ion- 
ized at least in part by the hard radiation field of the quasar, 
while the detection of the two carbon lines indicate that the 
gas is enriched in metals. As discussed by Humphrey et al. 
(2008b) and others, line ratios using only CIV A1549, Hell 
A1640 and CIII] A1908 do not provide strong constraints 
on the metallicity of the emitting gas: metallicities as low 
as ~5 per cent solar, or as high as a few times solar, are 
possible for the LAN. However, the non-detection of narrow 
NV emission argues against very high metallicities (see e.g. 
Vernet et al. 2001; Humphrey et al. 2008b). 

The Lya/Hell A1640 ratio, on the other hand, is a 



factor of ~2 higher than the average value, and is signifi- 
cantly higher than the range of values ~8-15 that 'normal' 
AGN-photoionization models can produce (Humphrey et al. 
2008b). Judging from the fitting carried out by van Ojik et 
al. (1997) for the extended Lya emission, the pre-absorbed 
Lya flux ought to be a factor of about xl.3 higher, which 
would raise the Lya/Hell ratio to ~29. Lya excesses of 
this kind have previously been identified in the LANs of 
a number of z>2 radio galaxies (e.g. van Ojik et al. 1994; 
Villar-Martm et al. 2007b), with the previously proposed ex- 
planations being (a) cooling radiation from low-metallicity 
infalling gas; (b) higher than expected nebular tempera- 
tures, due to low gas metallicity in the photoionized neb- 
ulae; (c) resonance flourescence; or (d) photoionization by 
young stars (see e.g. Villar-Martm et al. 2007b). In the case 
of the Lya nebula around TXS 1436+157, however, we must 
discard options (a) and (b) due to the substantial level of 
chemical enrichment of the nebula, so that only (c) and (d) 
remain as plausible explanations. Measuring the Lya/Ha 
ratio (e.g. Hayashi et al. 2012), or polarimetry of the Lya 
emission (e.g. Hayes, Scarlata & Siana 2011), would be use- 
ful to distinguish between these possibilities. 

4.7 The Nature of the Absorber and Lya nebula 

Assuming that the narrow Hell emission at the position of 
the quasar is emitted close to the systemic velocity, then 
the relative blueshift of the absorber implies it is in outflow. 
Various previous investigators have discussed the possible 
origins of similar blue-shifted, absorbing structures around 
high-z radio galaxies and quasars (e.g. van Ojik et al. 1997; 
Binette et al. 2000; Krause 2002; Jarvis et al. 2003; Wilman 
et al. 2004; Humphrey et al. 2008a). In the case of TXS 
1436+157, we favour a scenario adapted from the super- 
bubble scenario devised by Tenorio-Tagle et al. (1999). Dur- 
ing a powerful starburst, mechanical energy is deposited by 
supernovae explosions into the galaxian interstellar medium, 
producing in an expanding superbubble which decelerates as 
it sweeps up gas from the interstellar medium. The shell 
becomes Rayleigh- Taylor unstable and fragments (blows- 
out), thereby allowing the thermalized supernovae ejecta 
to rapidly escape into the low density halo of the galaxy, 
which then forms a new expanding shell sweeping up low 
metallicity gas from the halo. Left behind from the frag- 
mentation of the initial shell are Rayleigh- Taylor 'fingers' 
of cooler, denser, metal-enriched gas. These take ~100 Myr 
to fall from a radius of the order of ~10 kpc into the cen- 
tral regions of the galaxy, at a velocity of ~100 km s _1 (see 
S I4.6.1I) . to provide fuel to trigger the quasar and radio jets. 
The radiation field produced by the quasar then photoion- 
izes the infalling material, thereby allowing its detection in 
Lya and other lines, and also provides a background source 
for detection of the new shell, which has now expanded to 
a radius of several tens of kpc during the delay between the 
blow-out event and the triggering of the quasar. 

Interestingly, this scenario allows us to reconcile the 
simultaneous presence of infalling gas with roughly solar 
metallicity (Humphrey et al. 2007a, 2008b), which would 
correspond to the Rayleigh- Taylor 'fingers' left behind af- 
ter the blow-out of the initial shell, and an expanding shell 
of low-metallicity gas at larger radii (e.g. Binette et al. 
2000), corresponding to the second shell of swept-up mat- 
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ter. Naively, one may otherwise have expected gas falling 
into the galaxy to have substantially lower metallicity than 
gas flowing out. 

We must also consider whether it is energeticaly feasi- 
ble for the expansion of the absorbing shell to be starburst- 
driven. Averaged across the HI, CIV, NV and SilV lines (Ta- 
ble!]]), the current outflow velocity is ~140 km s -1 . Assum- 
ing this is a spherical shell of mass 

its kinetic energy would then be ^3.1xl0 58 erg. This serves 
as a lower limit on the mechanical energy deposited into the 
interstellar medium by supernovae. Using Fig. 116 of Lei- 
therer et al. (1999), we find that a modest star formation 
rate of ~20 Mq yr _1 , over a starburst lifetime of 100 Myr, 
would be sufficient to produce a mechanical energy luminos- 
ity of ~3xl0 58 erg. The total mass of stars formed in this 
starburst would be ~2 x 10 9 Mq , consistent with the stellar 
mass required to explain the alpha-element content of the 
shell 

Powerful radio sources can also provide a source of me- 
chanical energy to potentially drive an outflow (see e.g. Nes- 
vadba et al. 2006; Humphrey et al. 2010). Could the radio 
source have produced this outflow? The observed radius of 
the radio source is ~3"(Carilli et al. 1997). The fact that 
the radio source shows a relatively small arm-length asym- 
metry, shows a steep-spectrum and is lobe-dominated all 
imply that forshortening effects on its observed size can be 
neglected, so that its observed radius then corresponds to 
~25 kpc, which is substantially smaller than the rjs40 kpc 
of the absorber. While the ionizing radiation field of a quasar 
can also produce outflows (e.g. Humphrey et al. 2010), the 
kinetic energy of the shell around TXS 1436+157 is sub- 
stantially higher than even the most luminous quasars can 
reasonably inject (e.g. Feruglio et al. 2010). Thus, the ex- 
panding shell is unlikely to be driven by either the radio jets 
or the radiation field of the quasar. 

While the results presented in this paper directly con- 
cern only a single galaxy, they have wider implications for 
radio-loud active galaxies. Most radio-loud active galaxies 
at high-z have 10-100 kpc-scale Lya nebulae (e.g. Heckman 
et al. 1991; van Ojik et al. 1997) which are infalling where 
not perturbed by the radio jets (Humphrey et al. 2007) . A 
substantial fraction of these galaxies also show strong HI 
absorption features against the narrow Lya emission line, 
which are usually blueshifted relative to the emission (van 
Ojik et al. 1997), and whose spatial extension indicates a 
location outside of the extended Lya nebula. Clearly, then, 
TXS 1436+157 is not unique in this regard, and we argue 
that the scenario we have proposed above should also be ap- 
plicable to the wider population of high-z, radio-loud active 
galaxies. 

Finally, it is tempting to speculate that the shell 
may eventually escape from TXS 1436+157, and fragment 
into smaller self-gravitating baryonic entities, which at a 
later cosmic epoch may evolve into low-metallicity late-type 
galaxies. These galaxies would be chemically pre-enriched 
to a level 12 + log[0/H] ~7, in agreement with the low- 
metallicity theshold empirically found via observations of 
the most metal-poor star-forming galaxies currently known 
in the local Universe. Such galaxies would presumably be 
almost free from dark matter, somewhat analogous to high- 
metallicity 'tidal dwarf galaxies' ejected from galaxy mergers 
(see the review by Due 2010). 



5 CONCLUSIONS 

Using deep long-slit optical spectroscopy from the 10.4 
m GTC telescope, we have investigated the properties of 
narrow line emitting and absorbing gas associated with 
the z=2.54 radio-loud quasar TXS 1436+157. Our main 
conclusions are as follows. 

(i) The large Lya nebula around the quasar appears 
to be infalling onto the host galaxy at a rate of ~10-100 
Mq yr _1 . It shows anomalously strong Lya emission 
relative to Hell A1640, due to enhancement of Lya through 
resonance flourescence or through photoionization by stars 
(in addition to by the active nucleus). 

(ii) The absorber lies outside of the Lya nebula, and 
is likely to be a giant expanding shell enclosing the Lya 
nebula and the host galaxy. It has a hydrogen mass of 
> 1.6 x 10 11 M Q , a gas density of ^18 cm 3 , a geomet- 
rical thickness of ^18 pc, and a covering factor close to unity. 

(iii) The absorber is detected in CIV, NV and SilV 
lines, with low-ionization metal lines conspicuously absent, 
and is photoionized by the quasar. Using photoionization 
models we conclude that the metallicity of the absorbing 
gas is 12 + log[0/H] ^7.3 - much lower than that of the 
infalling Lya nebula. However, the detection of the metal 
lines shows that this is not pristine material. 

(iv) To explain the properties of the Lya nebula and 
the absorbing shell, we have proposed scenario in which a 
starburst-driven superbubble sweeps up material from the 
ISM; the bubble ruptures due to Ray leigh- Taylor instability, 
leaving behind filaments of cold, metal-enriched gas which 
fall back into the galaxy to fuel/trigger the quasar, while 
the hot supernova ejecta flows out into the low density halo 
and produces a new, larger shell of swept-up low metallicity 
gas. 

(v) We argue that when illuminated by a sufficiently 
luminous quasar, gaseous shells such as this should be 
detectable in emission lines. 
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